Quantum Error Correction for Quantum Memories 



B.M. Terhal 

Institute for Quantum Information 
RWTH Aachen University, 52056 Aachen, Germany 

^ February 15, 2013 

X> 

PL 

Abstract 

^— I This is a pedagogical review of the formahsm of stabihzer codes and their possible use in 

^ protecting quantum information. We discuss several codes of particular physical interest such 

i-C as the qubit-into-oscillator code, 2D topological error correction and 2D subsystem codes. The 

emphasis in this review is on the use of such codes as a quantum memory although we also 

discuss fault-tolerance theory and quantum universality. 
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1 Introduction 

When the idea of a quantum computer took hold in the '90s it was immediately realized that its 
implementation would require some form of robustness and error correction. Alexei Kitaev proposed 
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a scheme in which the physical representation of quantum information and reahzation of logical gates 
would be naturally robust due to the topological nature of the 2D physical system [1]. Around the 
same time Peter Shor formulated a first quantum error-correcting code and proved that a quantum 
computer could be made fault-tolerant [2] . Several authors then established the fault-tolerance thresh- 
old theorem (see Theorem [T]) which shows that in principle one can realize almost noise-free quantum 
computation using noisy components at the cost of moderate overhead. 

In this review we will discuss the basic ideas behind error correction: we will focus on quantum 
error correction using stabilizer codes for the purpose of making a quantum memory. Experimental 
efforts have not advanced very much into the domain of scalable quantum error correction. Several 
experiments exist of the 3-qubit (or 5-qubit) repetition code in liquid NMR, ion-trap, optical and 
superconducting qubits. Four qubit stabilizer pumping has been realized in ion-trap qubits . Some 
topological quantum error correction has been implemented with eight-photon cluster states in |1] 
and a continuous- variable version of Shor's 9-qubit code was implemented with optical beams [5]. 
Given the advances in coherence times and ideas of multi-qubit scalable design, in particular in ion- 
trap and superconducting qubits, one may hope to see more scalable error correction, fine-tuned to 
experimental capabilities and constraints, in the years to come. 

1.1 Shor's Code and The StabiUzer Formahsm 

The smallest classical code which can correct a single bit-flip error (represented by Pauli X is 
the 3-(qu)bit repetition code where we encode |o) = |000) and |l) = |111). A single error can be 
corrected by taking the majority of the three bit values and flipping the bit which is different from the 
majority. In quantum error correction we don't want to measure the 3 qubits to take a majority vote, 
as we would immediately loose the quantum information represented in the phase and amplitude of 
an encoded state j?/') = a |0) + /3 

But we can imagine measuring the parity checks ZxZ2 and ^'2.^3 without learning the state of each 
individual qubit. Fig. [T]^a) shows the quantum circuit which measures a parity check represented by 
a Pauli operator P using an ancilla qubit. Other than giving us parity information, the ideal parity 
measurement also provides a discretization of errors which is not naturally present in elementary 
quantum systems. Through interaction with classical or quantum systems the amplitude and phase 
of a qubit will fluctuate over time: bare quantum information encoded in atomic, photonic, spin 
or other single quantum systems is barely information as it is constantly at drift. An ideal parity 
measurement discretizes this continuum of errors into a discrete set of a Pauli errors (X, Y, Z or I on 
each qubit) which are amenable to correction. If the parity checks Z1Z2 and Z2Z3 have eigenvalues 
+1, one concludes no error. An outcome, say, of Z1Z2 = —1 and Z2Z3 = 1 is consistent with the 
erred state Xi I?/)) where is any encoded state. 

Let us informally introduce some of the notions used in describing a quantum (stabilizer) code. For 
a code C encoding k qubits, one defines k pairs of logical Pauli operators {Xi, Zi), i = 1, . . . k, such 
that XiZi = —ZiXi while logical Pauli operators with labels i and i' mutually commute (realizing 
the algebra of Pauli operators acting on k logical qubits). For the 3-qubit code we have X = X1X2X3 
and Z = Zi. 

The code space of a code C encoding k qubits is spanned by the codewords \x) where x is a fc- 
bitstring. All states in the codespace obey the parity checks, meaning that the parity check operators 
have eigenvalue +1 for all states in the code space. In other words, the parity checks act trivially 

^Pauli ax = X = ( ^ ^\a2 = Z= (^ -l) ^'^"^ '^^^"^^(i oj^ iXZ. 
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on the codcspacc. The logical operators of a quantum error-correcting code are non-unique as we 
can multiply them by the trivially-acting parity check operators to obtain equivalent operators. For 
example, Z for the 3-qubit code is either Zi or Z2, or Z3 or Z1Z2Z3. 

The 3-qubit repetition code docs not protect or detect Z (dephasing) errors as the parity checks 
only measure information in the Z-basis (Mz)- Shor's 9-qubit code was the first quantum error- 
correcting code which encodes a single qubit and corrects any single qubit Pauli error. Shor's code 
is obtained from the 3-qubit repetition code by concatenation. Code concatenation is a procedure 
in which we take the elementary qubits of the codewords of a code C and replace them by encoded 
qubits of a new code C . In Shor's construction we choose the first code C as the repetition code in 

the Hadamard-rotated basis {H : X ^ Z) with codewords |+) = |+ + +) and |— ) = | ). The 

parity checks of C are X1X2 and X2X3 and the logical operators are Zc = ZiZ2Z^ and Xc = Xi. 
As the second code C we choose the normal 3-qubit repetition code, i.e. we replace |-|-) by |+) = 
^(1000) + I111)) etc. 



(a) 



(b) 
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Figure 1: Measuring parity checks the quantum circuit-way. (a) Circuit to measure the ±1 eigenvalues 
of a unitary multi-qubit Pauli operator P. The gate is the controUed-P gate which applies P when 
the control qubit is 1 and / if the control qubit is 0. One measures the qubit in the |±) basis denoted 
as Mx- (b) Realizing the evolution exp{—i9P/2) itself (with Rx{0) = exp[—i6X/2). (c) Realization 
of circuit (a) using CNOTS when P — (d) Realization of circuit (a) using CNOTs when 

P = Z1Z2Z2Z4. 

We get all the parity checks for the concatenated 9-qubit code by taking all the parity checks 
of the codes C and taking the C'-encoded parity checks of C. For Shor's code this will give: the 
Z-checks Z1Z2, Z2Z3, Z4Z5, Z^Zq, Z^Zg and Z^Zg (from three uses of the code C) and the X-checks 
X4X5X6X7X8X9 (from the parity checks X1X2 and X2X3 where X is the logical 
operator of the code C). The non-unique logical operators of the encoded qubit are Z — Z1Z4Z7 and 
X = XiX2X^. 

This code can clearly correct any X error as it consists of three qubits each of which is encoded 
in the repetition code which can correct an X error. What happens if a single Z error occurs on any 
of the qubits? A single Z error will anti-commute with one of the parity X-checks or with both. For 
example, the error Zi anti-commutes with X1X2XZX4X5XQ so that the state Zi l'^) has eigenvalue 
— 1 with respect to this parity check. The eigenvalues of the parity check operators are called the 
error syndrome. Aside from detecting errors (finding —1 syndrome values) the error syndrome should 
allow one to infer which error occurred. 

One feature of quantum error correcting codes which is different from classical error correcting 
codes, is that this inference step does not necessarily has to point to a unique error. For example: 
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for the 9-qubit code, the error Z\ and the error Z2 have an equivalent effect on the codespace as 
Z\Z2 is a parity check which acts trivially on the code space. The syndromes for errors which are 
related by parity checks are always identical. The classical algorithm which processes the syndrome to 
infer an error, — this procedure is called decoding — , does not need to choose between such equivalent 
errors. But there is further ambiguity in the error syndrome. For Shor's code the error Z\ and the 
error Z^Z-j have an identical syndrome as Z\Za,Z-j is the Z operator which commutes with all parity 
checks. If we get a single non-trivial (—1) syndrome for the parity check X1X2X3X4X5X6 we could 
decide that the error is Z^ or But if we make a mistake in this decision and correct with 

while Z\ happened then we have effectively performed a Z without knowing it! This means that the 
decoding procedure should decide between errors, — all consistent with the error syndrome — , which 
are mutually related by logical operators. 

How is this decision made? We could assign a probability to each possible error: this assignment 
is captured by the error model. Then our decoding procedure can simply choose an error (or class 
of equivalent errors), consistent with the syndrome, which has highest probability. We will discuss 



the procedure of decoding more formally in Section 1.1.1 after we introduce the stabilizer formalism. 
For Shor's code we decode the syndrome by picking a single qubit error which is consistent with the 
syndrome. If a two-qubit error occurred we may thus have made a mistake. However, for Shor's code 
there are no two single-qubit errors E\ and E2 with the same syndrome whose product E\E2 is a logical 
operator. This implies that Shor's code can correct any single qubit error. It is a [[ri. A;, c/]] = [[9, 1, 3]] 
code, encoding k = \ qubit into n = 9 (n is called the block size of the code) and having distance 
d = ?). The distance d of the code is defined as the minimum weight of any logical operator (see the 
formal definition in Eq. ([T])). The weight of a Pauli operator is the number of qubits on which it acts 
non-trivially. We take the minimum weight as there are several logical operators and the weight of 
each one of them can be varied by multiplication with parity checks. It is simple to understand why 
a code with distance ci = 2t + 1 can correct t errors: errors of weight at most t have the property that 
their products have weight at most It < d and therefore the product can never be a logical operator 
as those have weight d or more. Thus our decoding procedure which picks an error of weight at most 
t can never lead to a logical error. If errors only takes place on some known subset of qubits, then a 
code with distance d can correct (errors on) subsets of size d—1 as the product of any two Pauli errors 
on this subset has weight at most d — 1. In other words, if c? — 1 or fewer qubits of the codeword fall 
into a black hole (or those qubits are lost or erased by other means) one can still recover the entire 
codeword from the remaining qubits. 

Clearly the usefulness of (quantum) error correction is directly related to the error model; it hinges 
on the assumption that low-weight errors are more likely than high-weight errors. Error-correcting 
a code which can perfectly correct errors with weight at most t, will lead to failure with probability 
roughly equal to the total probability of errors of weight larger than t. This probability for failure 
of error correction is called the logical error rate. The goal of quantum error correction is to use 
redundancy and correction to realize logical qubits with logical error rates below the error rate of the 
elementary constituent qubits. 

It seems rather simplistic to use error models which assign X, Z and Y errors probabilistically 
to qubits as in real quantum information amplitudes and phases are drifting rather than undergoing 
discretized errors. Ideal parity measurement can induce such discrete error model stated in terms of 
probabilities, but as parity measurements themselves will be inaccurate in a continuous fashion, such 
fully digitized picture is an oversimplification. The theory of quantum fault-tolerance, see Section 



1.3, has developed a framework which allows one to establish the results of quantum error correction 



and fault-tolerance for very general quantum dynamics obeying physical locality assumptions (see the 



4 



comprehensive results in |6j). However, for numerical studies of code performance it is impossible to 
simulate such more general open system dynamics and several simple error models are used to capture 
the expected performance of the codes. 

Three important remarks can be made with this general framework in mind. The simplest message 
is that a code which can correct t Pauli errors can in fact correct any possible error on a subset of t 
qubits (described for example by some noisy superoperator or master equation for the qubits). This 
point can be readily understood in the case of ideal parity check measurements where the measure- 
ment projects the set of possible errors onto the set of Pauli errors. Secondly, errors can be correlated 
in space and time arising from non-Markovian dynamics, but as long as (1) we use the proper estimate 
of the strength of the noise (which may involve using amplitudes and norms rather than probabili- 
ties) and (2) the noise is sufficiently short-ranged (meaning that noisy interactions between distant 
uncoupled qubits are sufficiently weak [7J), fault-tolerance threshold results can be established. The 
third remark is that qubit coding does not directly deal with leakage errors. As many elementary 
qubits are realized as two-level subspaces of higher- dimensional systems to which they can leak, other 
protective mechanisms such as cooling (or teleporting to a fresh qubit) will need to be employed in 
order to convert a leakage error into a regular error which can be corrected. 



Let us come back to Shor's code and imagine that the nine qubits are laid out in a 3 x 3 square 
array as in Fig. 13 It looks relatively simple to measure the parity Z-checks locally, while the weight-6 
X-checks would require a larger circuit. But why should there be such asymmetry between the X- 
and Z-checks? Imagine that instead of measuring X= i = X1X2X3X4X5X6, we measure (in parallel 
or sequentially) the eigenvalues of X1X4, X2X5 and X^Xq and take the product of these eigenvalues 
to obtain the eigenvalue of X= 1. The important property of these weight-2 operators is that all 
individually commute with the logical operators X and Z of the Shor code, hence measuring them 
does not change the expectation values of X and Z. These weight-2 X-checks do not commute with 
the weight-2 Z-checks however. If we first measure all the weight-2 X-checks and then measure 
the Z-checks, then with the second step the eigenvalues of individual X-checks are randomized but 
correlated. Namely, their product X1X2X3X4X5X6 remains fixed as X1X2X3X4X5X6 commutes with 
the weight-2 Z-checks. By symmetry the weight-2 X-checks commute with = Z1Z2Z4Z5Z7Z8 
and Z||^2 = ZiZ'^Z'^Z^Z^Z^. Viewing the Shor code this way we can imagine doing error correction 
and decoding using the stable commuting parity checks X=^i, X=^2, Z|| 1, Z|| 2 while we deduce their 
eigenvalues from measuring 12 weight-2 parity checks. Shor's code in this form is the smallest member 



in the family of Bacon-Shor codes [[n^, 1, n]] whose qubits can be laid out in a n x n array as in Fig. 13 



see Section 2.2 The Bacon-Shor code family in which non-commuting (low-weight) parity checks are 
measured in order to deduce the eigenvalues of commuting parity checks is an example of a (stabilizer) 
subsystem code. 



1.1.1 Stabilizer (Subsystem) Code^ 

Shor's code and many existing codes defined on qubits are examples of stabilizer codes [S]. Stabilizer 
codes are attractive as (i) they are the straightforward quantum generalization of classical binary 
linear codes, (ii) their logical operators and distance are easily determined, (iii) understanding of how 
to construct universal sets of logical gates and (iv) the numerical analysis of code performance are 
relatively simple. The main idea of stabilizer codes is to encode k logical qubits into n physical qubits 
using a subspace, the codespace, C C (C^)®*^ spanned by states \ip) that are invariant under the action 

^Readers less interested in this general framework can skip this section without major inconveniences. 
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of a stabilizer group S, 

Here S is an Abelian subgroup of the Pauli group Vn = Zi,. . . ,X„, Z„) such — / ^ S. For 

any stabihzer group S one can always choose a set of generators Si, ... , S^, i-e. S = {Si, . . . , Sm), 
such that Sa € Vn are hermitian Pauh operators. The mutually commuting parity checks which we 
considered before are the generators of the stabilizer group. If there are n — k linearly independent 
generators (parity checks) then the codespace C is 2'^-dimensional, or encodes k qubits. The weight 
\P\ of a Pauli operator P = Pi . . . P„ G "P^ is the number of non-identity single-qubit Pauli operators 
Pi. If the code encodes k logical qubits, it is always possible to find k pairs of logical operators 
{Xj, Zj)j=i^,,,^k- These logical operators commute with all the parity checks, i.e. they commute with 
all elements in S as they preserve the codespace. However they should not be generated by the parity 
checks themselves otherwise their action on the code space is trivial. The centralizer C{S) of S in Vn 
is defined as C{S) = {P G Pn|Vs G iS, Ps = sP}, i.e. all operators in Vn which commute with S. 
We thus have C{S) = {S, Xi, Zi, . . . , X^, Z^), i.e. the logical operators of the code are elements of 
C{S) \ as they are in C{S) but not in S. The distance ci of a stabilizer code can then be defined as 

d = min IPI. (1) 

P6C(5)\5 

Error correction proceeds by measuring the error syndrome s which are the ±1 eigenvalues of the 



generators of S. As we mentioned in Section 1.1 this syndrome will not point to a unique Pauli error 
but an equivalent class of errors. This class is a coset ES of the group S in Vn, where E = E{s) 
is an error which is consistent with the syndrome s. We can denote the coset as [E] = ES with 
coset representative Pj^and associate a total probability Prob([P]) = ^^g^. Prob(Ps) with the coset 
depending on some error model which assigns a probability Prob(P) to every Pauli operator P G Pn. 
If the syndrome is consistent with [E] it is also consistent with [EP] where P is any logical operator 
(as the logical operators commute with S). 

Maximum likelihood decoding is the procedure by which, given a syndrome s and a coset repre- 
sentative E{s), one finds the logical operator P (which could be /) which has the maximum value 
for Prob([PP]). It is important to consider how efficiently (in the number n of elementary qubits) 
this task can be done as Prob([PP]), — essentially a partition function — , is a sum over the number 
of elements in S which is exponential in n. For a simple depolarizing error model where each qubit 
undergoes a X, y or Z error with probability p/3 and no error with probability 1 — p, we have 
Prob([PP]) = (1 — p)" ^^g^ exp(— /3|PPs|) with inverse 'temperature' /3 = ln(3(l —p)/p). For small 
error rates p ^ 1 corresponding to low temperatures /3 — )• oo, the value of the partition function is 
dominated by the contribution of the lowest-weight error (s). Thus, instead of maximum likelihood 
decoding which compares the relative values of Prob([PP]), one can also opt for minimum- weight 
decoding. In minimum- weight decoding one simply picks an error P(s), consistent with the syndrome 
s, which has minimum weight |P|. We will discuss this decoding method for the surface code in 
Section H 

Subsystem stabilizer codes can be viewed as stabilizer codes in which some logical qubits, called 
gauge qubits, are not used to encode information [H] . The state of these qubits is irrelevant and can be 
fixed (gauge-fixing) or left variable. The presence of the gauge qubits sometimes lets one simplify the 
measurement of the stabilizer parity checks as the state of the gauge qubits is allowed to freely change 
under these measurements. One takes a stabilizer code S and splits its logical operators [X^, Zi) 
into two groups: the gauge qubit logical operators {Xi,Zi), i = l...m and the remaining logical 



^Note that left and right cosets are the sarae modulo trivial errors proportional to /. 
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operators (Xj, Zj) with i = m + 1, . . . k. We define a new subgroup Q = {S, Xi, Zi, . . . Xm, Zm) 
which is non-Abehan and contains S. As Q is non-Abehan we can consider its center, i.e. Q r\C{Q) = 
{P E Q\Wg & Pg = gP} = S (modulo trivial elements). 

If we measure the generators of the group Q we can deduce the eigenvalues of S. Since the k — m 
logical operators (Xj, Zi), i = m + 1, . . . ,k commute with Q, these logical operators are unaffected by 
the measurement. A priori there is no reason why measuring the generators of Q would be simpler than 
measuring the generators of the stabilizer S. In the interesting constructions such as the Bacon-Shor 
code and the subsystem surface code discussed in Section |2} we gain as we measure very low-weight 
parity checks in Q (while we lose by allowing more qubit-overhead or declining noise threshold). 

The distance of a subsystem code is not the same as that of a stabilizer code, Eq. ([T]), as we should 
only consider the minimum weight of the k — m logical operators. These logical operators are not 
unique as they can be multiplied by elements in S but also by the logical operators of the irrelevant 
gauge qubits. This motivates the definition of the distance as d = mmp^c{S)\g \P\- As errors on the 
gauge qubits are harmless it means that equivalent classes of errors are those related to each other 
by elements in Q. Given the eigenvalues of the stabilizer generators, the syndrome s, the decoding 
algorithm considers cosets E{s)Q in Vn, denoted as [E] = EQ. Maximum likelihood decoding proceeds 
by determining the coset [EP] which has a maximum value for Prob([-EP]) = Xlgee Prob(£'P(7) where 
P varies over the possible logical operators. 



1.1.2 Stabilizer Code Examples and The CSS Construction 

We discuss a few small examples of stabilizer codes to illustrate the formalism. For the two-qubit 
code with |o) = ^(|00) + |11)) and |T) = ^(|01) + |10)) we have X = Xi or X = X2 and Z = Z1Z2. 
The code can detect any single Z error as such error maps the two codewords onto the orthogonal 
states :^(|00) — |11)) and ;^(|01) — |10)) (as Z is of weight-2). The code can't detect single X errors 
as these are logical operators. 

The smallest non-trivial quantum code is the [[4,2,2]] error-detecting code. Its linearly indepen- 
dent parity checks are X1X2X3X4 and Z1Z2Z3Z4: the code encodes 4 — 2 = 2 qubits. You can verify 
that you can choose Xi = XiX2,Zi = Z1Z3 and X2 = X2X4, Z2 = Z3Z4 as the logical operators 
which commute with the parity checks. The code distance is 2 which means that the code cannot 
correct a single qubit error. The code can however still detect any single qubit error as any single 
qubit error anti-commutes with at least one of the parity checks which leads to a nontrivial —1 syn- 
drome. Alternatively, we can view this code as a subsystem code which has one logical qubit, say, 
qubit 1, and one gauge qubit, qubit 2. In that case Q = (X1X2X3X4, ZiZ2Z'3Z4, X2X4) = 
(Z'iZ'2, Z3Z4, X1X3, X2X4), showing that measuring weight-2 checks would suffice to detect single 
qubit errors on the encoded qubit 1. 

Another example is the stabilizer code Cq (defined in pU]) with parity checks X1X4X5X6, X1X2X3X6, 
Z1Z4Z5ZQ and Z1Z2Z3ZQ acting on 6 qubits. This code has 4 independent parity checks, hence it 
encodes 6 — 4 = 2 qubits with the logical operators Xi = X2X3, Zi = Z1Z2Z4 and X2 = X1X2X4, 
Z2 = -^4-2^5- As its distance is 2, it can only detect single X or Z errors (but note that it can correct 
a single X error on qubit 1, or Z error on qubit 2). 

One can concatenate this code Cq with the code [[4, 2, 2]] (called C4 in p^) by replacing the three 
pairs of qubits, i.e. the pairs (12), (34) and (56), in Cq by three sets of C4-encoded qubits, to obtain 
a new code. This code has thus n = 12 qubits and encodes k = 2 qubits. We can represent these 12 
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qubits as 3 sets of 4 qubits such that the X-checks read 
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and the logical operators are 
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One can verify that the minimum weight of the logical operators of this concatenated code is 4. Thus 
the code is a [[12, 2,4]] code, able to correct any single error and to detect any three errors. 



One could repeat the concatenation step and recursively concatenate Cg with itself (replacing a 
pair of qubits by three pairs of qubits etc.) as in Knill's C^/Cq architecture [TU] or, alternatively, 
recursively concatenate C4 with itself as was considered in [TT]. Note that in general when we con- 
catenate a [[ni,l,(ii]] code with a [[n2, 1,(^2]] code, we obtain a code which encodes one qubit into 
n = nin2 qubits and has distance d = did2- Code concatenation is a useful way to obtain a large 
code from smaller codes as the number of syndrome collections scales linearly with the number of 
concatenation steps while the number of qubits and the distance grows exponentially with the number 
of concatenation steps. In addition, decoding of a concatenated code is efficient in the block-size n 
of the code and the performance of decoding can be strongly enhanced by using message passing 
between concatenation layers |12j . 

The smallest stabilizer code which encodes 1 qubit and corrects 1 error is the [[5,1,3]] code; you 
can find its parity checks in [13]. Another well-known code is Steane's 7-qubit code, [[7, 1,3]] which is 
constructed from two classical codes using the Calderbank-Shor-Steane (CSS) construction. Classical 
binary linear codes are fully characterized by their parity check matrix H. The parity check matrix 
Hi of a code Ci encoding ki bits is a (n — /ci) x n matrix with 0,1 entries where linearly independent 
rows represent the parity checks. The binary vectors c G {0,1}" which obey the parity checks, i.e. 
He = (where addition is modulo 2), are the codewords. The distance ci = 2t + 1 of such classical 
code is the minimum (Hamming) weight of any codeword and the code can correct t errors. 

We can represent a row r of Hi of a code Ci by a parity check operator s{Z) such that for the 
bit Tj = 1 we take s{Z)i = Z and for bit rj = we set s{Z)i = I. These parity checks generate some 
stabilizer group Si{Z). In order to make this into a quantum code with distance larger than one, one 
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needs to add X-type parity checks. These could simply be obtained from the {n — ^2) x n parity check 
matrix H2 of another classical code C2. We obtain the stabilizer parity checks S2{X) by replacing the 
Is in each row of this matrix by Pauli X and / otherwise. But in order for S = {Si{Z) , S2{X)) to 
be an Abelian group the checks all have to commute. This implies that every parity X-check should 
overlap on an even number of qubits with every parity Z-check. In coding words it means that the 
rows of H2 have to be orthogonal to the rows of Hi. This in turn can be expressed as C Ci where 
C2 is the code dual to C2 (codewords of are all the binary vectors orthogonal to all codewords 
CGC2). 

In total S = {Si{Z) , S2{X)) will be generated by 2n — ki — k2 independent parity checks so that 
the quantum code encodes ki + k2 — n qubits. The distance of the quantum code is the minimum of 
the distance d{Ci) and d{C2) as one code is used to correct Z errors and the other code is used to 
correct X errors. 

A good example is Steane's code which is constructed using a classical binary code C which 
encodes 4 bits into 7 bits and has distance 3. Its parity check matrix is 

/ 1 1 1 1 \ 
= 1 1 1 1 . (2) 
\ 1 1 1 1 / 

The codewords c which obey He = are linear combinations of the 7 — 3 = 4 binary vectors 
(1, 1, 1, 0, 0, 0, 0), (0, 0, 0, 1, 1, 1, 1), (0, 1, 1, 0, 0, 1, 1), (1, 0, 1, 0, 1, 0, 1) where the last three are the rows 
of the parity check matrix: these are also codewords of C"*-. Hence C"*- C C and we can use the CSS 
construction with Ci = C and C2 = C to get a quantum code. As C*-*" (as well as C) has distance 
3, the quantum code will have distance 3 and encodes one qubit. The parity checks are Z^Z^Z^Z^, 
Z2Z3ZQZJ, ZiZ^Z^Zj and X^X^XqXj, X2X3X6X7, X1X3X5X7. 



1.2 QEC Conditions and Other Small Codes of Physical Interest 

One may ask what properties a general (not necessarily stabilizer) quantum code, — defined as some 
subspace C of a physical state space — , should obey in order for a certain set of errors to be correctable. 
These properties are expressed as the quantum error-correcting conditions which can hold exactly 
or only approximately. We encode some k qubits into a codespace C so that |z) are the codewords 
encoding the fc-bit strings i. Assume there is a set of errors S = {Ei} against which we wish to correct. 
For example we can describe the short-term decoherence effects on the qubits by a superoperator 
S{p) = YliEipEj which has the Kraus error operators Ei. The quantum error-correcting conditions 
say that there exists an error-correcting operation, a reversal of the error, if and only if the following 
conditions are obeyed for all errors Ek,Ei E £ 

,{J\ElEi\j) = CMSij. (3) 

Here Cki is a constant independent of the codeword The condition for i = j says that the codewords 
are not distinguished by the error observables. The condition for i j indicates that the orthogonal 
codewords need to remain orthogonal after the action of the errors (otherwise we could not undo the 
effect of the noise). Note that if a code can correct the error set {Ei} it can also correct an error 
set {Fj} where each Fj is any linear combination of the elements Ei as the set {Fj} will also obey 
the quantum error-correcting conditions of Eq. ([3]). These QEC conditions can be generalized to the 
unified framework of operator quantum error correction [HI [15] which covers both subsystem codes as 
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well as error-avoidance techniques via the use of decoherence-free subspaces and noise-free subsystems. 



Quite commonly one describes an open system dynamics by a Markovian master equation of 
Lindblad form 

f^=-i[H,p]+C{p)^CUp) (4) 

where £(p) = J2j-^jP-^] " \{^\^j^P} "with quantum-jump or Lindblad operators Lj while if is a 
harmless or desired driving Hamiltonian. For short times r we have p(r) = S-r^pifi)) = EqpEq + 
J^iJ^iP-^i "with Eq ^ I — irH — ^T^^L\Li = I — 0(r) and Ei ^ y/rLi. Thus the error set is 
approximately given by the quantum jump operators Lj. 

Let us consider two simple codes which approximately obey the conditions in Eq. ([s]) for the phys- 
ically relevant decoherence due to (amplitude) damping for qubits and bosonic modes. Even though 
the [[5, 1, 3]] code is the smallest code which can correct against any single qubit error one can use 4 
qubits to approximately correct any amplitude-damping error which can model energy loss [16j. The 
noise process for amplitude damping on a single qubit is given by the superoperator S{p) = AipA\ 
with Ao = + y/T^)I + (1 - VI - k)Z) ^ I - 0{k) and Ai = y/na^. The codewords for the 
4-qubit amplitude damping code are |0) = 4^(|0000) + |1111)) and |T) = ^(|0011) + |1100)). Each 
qubit in this code is subjected to amplitude-damping noise. We wish to approximately correct against 
the error set Eo = Ei= Ai(g) A®^ Eg = ® O A^^, E^ = A®^ ® O and E^ = ® Ai 
which correspond to no damping and single qubit damping on any of the four qubits respectively. The 
authors in [IB] show that this code obeys the QEC conditions approximately with O(k^) corrections 
which is a quadratic improvement over the basic error rate k. 

Another interesting example is that of a single bosonic mode (with creation and annihilation 
operators a\a) which is used to encode a qubit in two orthogonal states which are approximately 
protected against photon loss [IT] . The damping process can be modeled with the Lindblad equation, 
Eq. (jij), with L = y/i^a while H = (u + |)a^a (which we can transform away by going to the rotating 
frame at frequency u). One can choose two Schrodinger cat states as encoded states 

|o+> = ^d") + 1-"))' |T+> = ^(R + H«))- (5) 



Here |a) is a coherent state |a) = exp(— |ap/2) |n) and N± = 2(1 ± exp(— 2|ap)) ^ 2. 

For sufficiently large photon number (n) = the states \±a) ,\±ia) (and thus |0+) and |l+)) 

are approximately orthogonal. The creation and manipulation of such cat states has been actively 
explored for cavity modes in micro-wave cavities, see e.g. [IH]. The code states are chosen such 
that loss of a photon from the cavity maps the states onto (approximately) orthogonal states. As 
a\a) = a \a), we have 

a|0+) = av^i\^Ji\L |0_), a|T+) =mv^iV^7^|T-> , (6) 
with |0_) = I — {\a) — |-a)) and |l_) = —T={\ia) — \-ia)). As we know the preservation of orthog- 



onality is a prerequisite for these code states to be correctable. One can verify that in the limit of 



\a\ 



— )■ oo we obey the QEC conditions^ Eq. dsl), for Eq = y/na and Ei = I — ^a'^a. The cat states 



"'if we were to use two coherent states as code states, say, |0) = \a) and |l) = |-a), the QEC conditions would not 
be obeyed, as {a I eIEo \a) ^ (-a| ElE^i \-a) for any a. 
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will gradually loose their intensity and thus their approximate protection under these processes as 
the vacuum state |0) is the only stationary state of the Lindblad equation. Counterintuitively they 
do this by the process with Kraus operator Ei, i.e. by not observing the photon emission, while the 
Kraus operator Eq maps the state |0+) onto |0_) and vice- versa |18]. The codespace (spanned by 
|0+) , is distinguished from the orthogonal erred space (spanned by |0_) and |l-)) by the parity 
operator exp^ina'^a) = J2ni~^)"' \^) ~ Peven — Podd- This parity operator has +1 eigenvalue for 
the states |0+) , and —1 eigenvalue for the states |0_) , By continuously monitoring the 

value of the parity operator exp(i7ra"''a) one could track the occurrence of errors. Better even would 
be the realization of a restoring operation which puts back an erred state with decayed amplitude 
Q,g-Kt/2 jj^i^Q code space while restoring the amplitude back to a as proposed in [T7] In [T7] the 
interaction with superconducting qubits coupled to the cavity is proposed to be used for encoding, 
correction and decoding of the cat states. 

One can generalize the stabilizer formalism to continuous-variable systems characterized by an 
infinite-dimensional Hilbert space [121 120] (see also [21] )• Of particular interest are codes which 
encode a discrete amount of information, a qubit say, in a harmonic oscillator [22] • Given are two 
conjugate variables p and q which represent generalized momentum and position, obeying = i. 
The idea is to encode the information such that small shifts in position or momentum correspond to 
correctable errors while logical operators are represented as large shifts. For a harmonic oscillator 
space, the Pauli group P„ can be generahzed to the Weyl-Heisenberg group generated by the unitary 
operators exp(itp) and exp(zsg) for real s and t. In order to define a qubit in this infinite-dimensional 
space we select a set of commuting check generators whose -|-1 eigenvalue space is two-dimensional. 

We will consider three examples. In our first example the codespace is a single state and essentially 
represents a quantum rotor. We choose Sq = (?^^ and S.p = e~^'^^ as commuting check operators |^ 
When Sp = 1 and Sq = 1 the eigenvalues of p are even integers while q should have eigenvalues 
mod TT. Defining n = p/2 and = 2q, we obtain a quantum rotor with conjugate variables 
n = 0, ±1, . . . and 27r-periodic phase 0. A realization of the quantum rotor is the quantization of a 
superconducting circuit where is the superconducting phase (across a junction) and h represents 
the number of Cooper pairs. As there is a unique state with fixed (say +1) eigenvalues for Sq and 
Sp, superconducting qubits such as the transmon qubit [23] use superpositions of states with different 
eigenvalues for 5*^ = e*"^ in order to encode information. 

A different choice of Sq and Sp leads to a real code which encodes a single qubit and has built-in 
protection. We choose as checks the operators Sq = e^^^ and Sp = e~^*'^^. Fixing the eigenvalues 
of these operators to be +1 leads to the discretization p = 0,±1,±2... and again q should have 
eigenvalues which are multiples of tt. Now there are two operators which commute with Sq and 
Sp but which mutually anti-commute: these are Z = e^^ and X = e~™^. One can verify the proper 
commutation relations of Sq, Sp, X and Z , using e^^e^ = e^^'^^e^e^ (when A, B are linear combinations 
of q and p). 

The state |o) (defined by Z |o) = |o) and Sp |o) = |o)) is a uniform superposition of states with 
g = 0, ±27r, . . .. Similarly, |l) corresponds to a uniform superposition of g = ±7r, iSvr, . . ., see Fig. |2] 
with a = vr. Consider the effect of shifts of the form e^^'^ where \5\ < 7r/2, which are correctable. Such 

^Such restorative process is physically impossible for all states in the code space as it would increase the distin- 
guishability of states by decreasing their overlap. 

^As q and p are dimensionfull quantities one should take Sq = e^^^l^° and — e^^^vlva where and set the 
scale so that g/gp etc. are dimensionless. In what follows we tacitly assume the existence of such a scale so that if, say, 
the codespace corresponds to integer f> it is an integer in units of po- 
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shifts map the codewords outside of the codespace as they do not commute with the stabihzer operator 
Sq. Error correction thus takes place by measuring q mod tt and applying the smallest shift which 
resets q = mod n. Similarly, the |+) is a uniform superposition of states with p = 0, ±2, ±4, . . . 
while |— ) is a uniform superposition of states with p = ±1, ±3, . . ., see Fig. [2} The qubit is protected 
against shifts e*'^'^ with |e| < 1/2. 

This code is realized as a Majorana fermion qubit where p = n counts the total number of electrons 
while g = is the vr-periodic conjugate phase variable. The |+) eigenstate of X with an even number 
of electrons correspond to the Majorana mode unoccupied while |— ) is the state with an odd number 
of electrons as the Majorana mode is occupied. The protection of the Majorana fermion qubit can 
thus also be understood from this coding perspective although the formalism does of course not shed 
light on how to realize this qubit. Another realization of this code, which does not use Majorana 
fermion qubits, but superconducting circuits is the O-vr qubit (see e.g. [21]) which is designed such 
that the superconducting phase difference between terminals has degenerate energy minima at and 
IT corresponding to the approximate codewords |o) and |l). 
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Figure 2: Picture from [22]: Amplitude of codewords for the stabilizer code with commuting checks 
Sq{a) = e^™^!'^ and S'p(a) = e~^*^" which encodes a qubit in an oscillator. 

More generally, we can parametrize this code by a real number a by taking the stabilizer checks 
as Sq = e^*'^'^/" and Sp = e~^*^" (above we took a = tc). The logical operators are Z = e^'^^l'^ and 
X = e~*^" [22], see the codewords in Fig. [2j The code can correct against shifts e*'^'' with |e| < ^ 
and e~*^^ where \S\ < |. One can use this code for encoding a qubit in a bosonic mode where q and 
p arise as quadrature variables, i.e. q = -^{a'^ + a) and p = -^{a^ — a) . 

The free Hamiltonian Hq = u^a^a + |) will periodically transform q into p and vice versa so it is 
natural to let Sg be of the same form as Sp and choose a = a/vt. The codewords of this code in Fig. [2] 
are not physically achievable as it requires an infinite amount of squeezing to prepare (superpositions 
of) of a quadrature eigenstates such as \q) or \p). Ref. [22] proposed to use approximate codewords: 
for example the approximate codeword |o) is a superposition of Gaussian peaks in g-space, — each 
one centered at integer multiples of 2^/7^ with width A — , in a total Gaussian envelope of width 1/k. 
Viewed as a superposition of p-eigenstates, such state is a superposition of peaks with width k and 
total envelope of width A~^. An error analyis of this approximate encoding was done in [23], while 
[2n] considered the preparation of the encoded states using cat states as in Eq. (|5]), squeezing and 



homodyne detection. In Section 2.1.3 we will consider a version of the surface or toric code which 



encodes an oscillator in a 2D coupled array of harmonic oscillators. 
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There are of course many codes which are not captured by the qubit (or qudit) stabihzer formahsm. 
Here I would hke to briefly mention the class of 2D topological qubit codes where the stabilizer 
checks are still commuting, but they are no longer simple Pauli operators. As Hamiltonians these 
correspond to the so-called 2D Levin- Wen models [2^, as codes they are called Turaev-Viro codes 
[28] . The advantage of these codes which generalize the 2D surface code in Section |2| is that universal 
quantum computation can achieved by purely topological means. The disadvantage from the coding 
perspective is that (1) the stabilizer checks are more complicated as operators, e.g. for the so-called 
Fibonacci code on a hexagonal lattice, the stabilizer checks act on 3 and 12 qubits and (2) no one has 
established a scheme for error correction or a noise-threshold for these codes. 



1.3 Error Correction and Fault Tolerance 

We have understood from the previous sections that the crucial element of quantum error correction 
for stabilizer codes is the realization of the (parity) check measurement as in Fig. [l] The immediate 
problem is that the parity check measurement suffers from the same imperfections and noise as any 
other gate or measurement that one may wish to do. The parity check measurement should be a 
non-demolition, strong single-shot, measurement which has the effect of projecting the state in either 
the code space or a discrete set of error spaces while giving at the same time accurate information 
about which event occurred. In practice a parity check measurement may arise as a continuous weak 
measurement leaving a classical stochastic data record which hovers around the value +1 (pointing to 
the state being in the code space) while occasionally jumping to a value centered around —1. One can 
imagine that this continuously acquired record is immediately fed back to unitarily steer the qubits to 
the code space |29]. Advantages of such immediate feed-back or 'closed-loop' form of error correction 
are that it eliminates the need for a fast and accurate quantum measurement in which weak quantum 
signals are amplified to our classical macroscopic world. A disadvantage which we will discuss at the 



end of Section 2A_ is that the correction mechanism is purely local in time and space. 

We could set up such local error correction mechanism by explicitly engineering a dissipative 
dynamics which drives/corrects the qubits towards the code space as proposed in [3l [30]. We assume 
that the open-system dynamics of code qubits and environment is described by a Lindblad equation 
as in Eq. Q. For simplicity we consider the case in which we would like to pump or drive four 
qubits into a state with even parity so that the 4-qubit parity Z-check, Z1Z2Z3Z4 has eigenvalue 
-fl. Imagine that we can engineer the dissipation (in the interaction picture) such that there is a 
single quantum jump operator L = -^XiPom with Pom = \{I — ^1^2^3^4)5 the projector onto 
the odd parity space, and H oc ZxZ2Z-iZ/^. Integration of the Lindblad equation gives rise to the 
time-dynamics p(t) = exp(t£tot)(p(t = 0)) with stationary states p determined by Ctot{.p) = 0. States 
supported on the even-parity subspace are dark states with £(p) = and [-ff, p] = 0. The odd-parity 
subspace is not stationary as the quantum jump operator L flips the first qubit so that an odd parity 
state becomes an even parity state pumping the system towards the stationary dark subspace. 

In pD] one considers the following stroboscopic evolution using an ancillary dissipative qubit or 
mode which approximately gives rise to such Lindblad equation. The idea is to alternate (or trotterize) 
the coherent evolution with H and the dissipative evolution with C for short periods of time r so that 
exp(r£iot) ^ exp(— zr[i7, .]) exp(r£). The dynamics of H can be obtained by a small modification of 
the parity check measurement circuits in Fig. [T| for the evolution exp(— i6'P) where P is a multi-qubit 
Pauli operator we can use the circuit in Fig. ^ib). 

The dissipative evolution C could be implemented for short times r <^ 1 using a circuit consisting 
of a dissipative ancilla coupled to the four qubits as in Fig. ^d). Instead of immediately measuring 
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the ancilla qubit, we apply a CNOT with the ancilla qubit as control and qubit 1 as target (to change 
the parity of the odd states to even), followed by decay (cr_) of the ancilla qubit from |1) to |0). These 
ideas of stabilizer pumping were experimentally tested on two and four ion-trap qubits in Ref. [5]. 



G 



Mr 

Figure 3: Concatenation: a gate G in the circuit M^-i is replaced by a rectangle consisting of error- 
correcting steps {E) followed by a fault-tolerant encoded realization of the gate. The process can be 
repeated for every elementary gate in the new circuit Mr. 

As any realization, closed or open-loop, of quantum error correction will suffer from inaccuracies 
there is no guarantee that one will improve coherence times by encoding a qubit in a code as it may 
introduce more errors that it takes away. This has been the reason for the construction of a theory of 
quantum fault tolerance which addresses this problem using recursively applied code concatenation. 
For simplicity we assume that every such elementary gate, idling step or measurement, — these are 
called locations in the circuit — , can fail independently with some error probability p (independent 
stochastic noise). In a concatenation scheme every qubit and operation in a quantum circuit is 
replaced by an encoded qubit resp. encoded operation and the process is recursively repeated. The 
encoded operation consists of an error correction step and a fault-tolerant realization of the operation, 
see Fig. [3| which together constitute a rectangle. For a code such as Steane's [[7, 1, 3]] code which can 
correct a single error, the fault-tolerance of the rectangle should be such that a single error in any 
of the locations of the rectangle cannot lead to two (incorrectable) errors in one code block. Then, 
if the elementary error rate scales as p, it follows that the encoded error rate scales as Cp^ as two 
elementary errors are required for a logical error. Here C is a constant which roughly counts the 
number of pairs of locations in the rectangle where failure can lead to a logical error. If Cp^ < p the 
concatenation step helps and r steps of concatenation will drive down the error rate to ~ p^*^ while 
the overhead in terms of qubits and gates increase only exponentially in r. The equality Cp^ = p sets 
the noise threshold Pc. 

This idea of repeated code concatenation was used in the early days of quantum error correction 
to prove the Threshold Theorem which says that fault-tolerant computation is possible with arbitrary 
small error rate if one is willing to tolerate an overhead which scales polylogarithmically with the size 
N of the computation to be performed (the size of a quantum circuit is the number of locations in 
it). 

Theorem 1 ([Si 1311 132L I33j) An ideal circuit of size N can be simulated with arbitrary small error 
S by a noisy quantum circuit subjected to independent stochastic noise of strength less than p < Pc 
where the noisy quantum circuit is of size 0{N{log Ny) with some constant c. 

The important question is of course what is the noise threshold pc and how large is the constant c? 
The best performing concatenated coding scheme to date is the C^/Cq scheme of Knill [10] for which 
he estimated a noise threshold as high as Pc ~ 3% but at the cost of huge overheads. In general. 
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the overhead that one incurs by code concatenation tends to be worse than the overhead (or mere 
complexity) that is incurred using the idea of topological error correction (Section [2]). One reason 
is that in topological quantum error correction parity check measurements are simply made robust 
by repeating the measurement needing no additional qubits. The constant c in the Theorem equals 
c = log2 S where S is the number of locations in a rectangle. For example, in a 2D realization of 
the concatenated Steane [[7, 1, 3]] code, the fault-tolerant CNOT has 5* = O(IO^) J3^ so that c ~ 10 
demonstrating the potential inefficiency of code concatenation. 

1.4 Universality 

An important feature of error correction with stabilizer (subsystem) codes is the fact that one never 
physically has to do the correcting Pauli operation as long as the qubits only undergo Clifford group op- 
erations. The decoding procedure gives a Pauli error which is interpreted as a frame, the so-called Pauli 
frame which we can efficiently track during the computation using the Knill-Gottesman theorem. 
The Clifford group C„ is a finite subgroup of the group of unitary transformations W(2") on n qubits. 
It is defined as the normalizer of the Pauh group: Cn = {U e W(2")|VP E P„, 3P', UPW = P'}. The 
generators of the Clifford group are the 2-qubit CNOT gate, the Hadamard H gate, the phase gate S 
and the Pauli's X, Z. The Knill-Gottesman theorem [33] proves that one can efficiently classically 
simulate any quantum circuit with Clifford group gates. One does this by tracking the stabilizer 
group, more precisely its generators, which has the input state of the quantum circuit as its unique 
+1 state. Every Clifford gate and measurement maps the stabilizer generators onto new stabilizer 
generators providing an efficient representation of the action of the quantum circuit. The presence of 
additional Pauli errors can thus be represented as additional updates of the stabilizer generators. 

For universal quantum computation one needs additional gates such as the T gate (tt/S rotation). 
Examples of universal gate sets are {if, T, CNOT}, {if, Toffoli} and {H,A{S)} where A{S) is the 
two-qubit controUed-S* gate |^ Even though Clifford group gates have no quantum computational 
power they can be used to develop a quantum substrate on which to build universal computation 
using stabilizer codes. This comes about by combining the following sets of ideas. First of all, 
note that stabilizer error correction by itself only uses CNOT gates, preparations of |+) , |— ) , |0) , |1) 
and measurements in the Z- and X-basis as is clear from Fig. [l| The T, A(S) and the Toffoli 
gate, each of which can be used with Clifford gates to get universality, are special unitary gates as 
they map Pauli errors onto elements of the Clifford group. One can define a Clifford hierarchy [36] 
C{j) = {U e W(2")|?7P„f/t C C{j - 1)} such that C(l) = P„, C(2) = C„. The T,A{S) and Toffoh are 
thus members of C(3). Such gates in C(3) (and similarly gates in C{j) for j > 3) can be realized with 
ancillas and Clifford group gates using quantum teleportation ideas [351 EZ] • The idea is illustrated in 
Fig. [5] for the T gate. One teleports the qubit on which the T gate has to act, prior to applying the 
gate, using the bottom one-bit teleportation circuit in Fig. |4j Inserting TT"^ prior to the corrective 
Pauli X and using that TXTt = e-'^'/^SX which is an element of the Clifford group, we obtain the 
circuit in Fig. [5j 

Example: We can do the same trick for the S gate, that is, we can reduce the S gate to the 
preparation of a \+i) = :^(|0) + i |1)) ancilla, a CNOT gate and a corrective Pauli Y. We get this 
from starting with the bottom circuit in Fig. |4]to which we apply the S gate at the output. We insert 

'"-Ml ">^^(i 

8 A{S) |6i, 62) = 5"^ I62) for 61, 62 = 0, 1. 
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Figure 4: Picture from [37] : so-called one-bit teleportation circuits. The measurement denoted by the 
meter is a measurement in the Z-basis and determines whether to do a Pauli on the output qubit. 
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Figure 5: Picture from [37] : using the ancilla T |+) in the dashed box, one can realize the T gate by 
doing a corrective operation SX. 



SS'' before the corrective Pauli X and use that SXS'' oc Y. The gate 5* can be commuted through 
the control-line of the CNOT as both gates are diagonal in the Z-basis on the control qubit. We thus 
need the ancilla SH |0) = ^(|0) + i |1)). 

How do we realize a universal set of logical fault-tolerant gates for a code? In principle, fault- 
tolerant gate constructions can be made for any stabilizer code [S]. The question is how to do 
computation with resource requirements and overheads similar as those needed for a quantum memory. 
Ideally, the computation threshold, i.e. the performance of the code when used for computation is 
close to the memory noise threshold, the performance of the code as a pure quantum memory. An 
example of a gate which does not require new resources is a transversal CNOT between two code 
blocks in which every qubit in the block is paired with a qubit in the other block in a CNOT gate. 
Such transversal CNOT represents a valid logical CNOT for stabilizer codes S = {S{X),S{Z)), i.e. 
the parity X-checks are the same as parity Z-checks. Examples are Steane's [[7, 1,3]] code or any 
other CSS code constructed from a classical code C with C"*- C C. 

In pS] it was shown for any 2D stabilizer code that the logical gates which can be performed by 
constant-depth circuits employing only local gates, are members of the Clifford group. The reason to 
focus on constant-depth local circuits is that such circuits are small and are naturally fault-tolerant. 
Namely, any number of errors which occurs in the circuit will only affect a patch of 0(1) qubits on the 
2D lattice, and such 0(1) error patches are correctable. Hence we expect that their implementation 
does not negatively impact the noise threshold or the overhead. 

As we will discuss in Section [2] one can implement the logical H and CNOT by topological means 
in the 2D surface code, but one cannot do the logical S gate in this manner. For the surface code one 
can do the logical S in the same fashion as the logical T, see below (for a different logical S trick, see 
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[39]). Other stabilizer codes, 2D color codes, have been found which allow for a topological realization 
of the full Clifford group |10j. Interestingly, Ref. [38] shows that for D-dimensional stabilizer codes 
all logical encoded gates which are composed from constant-depth logical gates are contained in the 
Clifford hierarchy C{D). This is a subtle result as we can realize a fault-tolerant set of universal gates 
for any stabilizer code, but apparently we cannot do this by composing a sequence of constant-depth 
encoded gates. 

For gates such as the T gate the method of magic-state-distillation has been developed jH] show- 
ing how to realize these gates fault-tolerantly assuming the availability of noiseless Clifford group 
operations. Thus, once we have built a low-noise Clifford computation substrate, universal quantum 
computation can be bootstrapped from it. The strength of the distillation scheme is that the noise 
threshold for it to work is high, approximately 15%. The downside is that the qubit/gate overhead 
per logical T gate is orders of magnitude larger that of a 'topological' CNOT (see e.g. Fig. 11 in 
|12]). In a nutshell, the ideas are as follows. We implement the T gate at the logical level using Fig. |5] 
which requires the preparation of low- noise logical ancillas T\+). We can obtain such an ancilla 
in non-fault-tolerant noisy manner by, for example, injecting several noisy unencoded ancillas into 
the code [10]. From many of these noisy encoded ancillas we distill using logical H, CNOTs and 
measurements, a single low-error encoded ancilla. 

Note that prior to the action of logical non-Clifford gates one does need to physically implement 
the inferred logical Pauli errors obtained from decoding as gates such as the T gate do not map Pauli 
errors onto Pauli errors but onto Clifford errors. What one never needs to physically do is a correction 
which maps one back to the +1 eigenspace of the stabilizer S: any syndrome eigenspace of 5 is a 
good code. 

The upshot of these considerations is that 2D and 3D stabilizer codes will be most suitable for 
building a quantum memory and performing Clifford group operations. The goal of universal quantum 
computation within the same platform can be reached using methods such as injection-and-distillation 
but the additional overhead and complexity of distillation (and demands for fast decoding) are consid- 
erable. Several recent papers have been devoted to reducing this overhead, see e.g. [13] and references 
therein. 



2 2D (Topological) Error Correction 

In this section we discuss three stabilizer/subsystem codes in which the parity checks act locally on 
qubits laid out on a 2D lattice. Two of these codes are topological codes: if the code was defined on 
a 2D lattice with periodic boundary conditions, — a torus — , then the logical operators would relate 
to the non-trivial loops around the torus. 

Most experimental qubit quality is assessed by the twin parameters T2 and Ti which capture the 
decay of off-diagonal matrix elements resp. the relaxation to steady state of the diagonal matrix 
elements in the qubit density matrix. A simple estimate for an error rate can be made by taking 
r/Ti to be the probability for an X error where r is a time-scale which captures the time needed for 
error-correction through single and two-qubit gates. Similarly T/T2 can provide an estimate for the 
Z error rate. In numerical or analytical studies of code performance, one uses simple error models 
such the independent depolarizing noise model. Independent depolarizing noise assumes that every 
qubit independently undergoes a X,Y or Z error with equal probabilities p/3 and no error with 

Z?-dimensional stabilizer code is a stabilizer code for which the parity checks are local on some D-dimensional 

lattice. 
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probability 1 — p. A related noise model is that of independent X and Z errors in which a qubit 
can undergo independently an X error with probability p and a Z error with probability p in each 
time-step. In all three codes the parity checks are either X or Z-like, detecting either Z or X errors. 
In addition, the parity Z- and X-checks have the same form so error correction proceeds identically 
(but simultaneously) for both types of errors. 

We will consider codes which encode a single qubit in a block of n qubits with n = 0{L^) with 
L the linear size of the 2D array. Important parameters of the code performance are the pseudo- 
threshold Pc{L) which is the error probability at which the elementary error rate p equals the logical 
error rate for a given size L, the asymptotic threshold Pc = lim„_^ooPc(-^) and the logical error rate 
p{p, L) which is a function of the block size and the elementary error rate p. For the Bacon-Shor 
code in Section 2.2, the asymptotic threshold p^ = 0, hence it is of interest to consider what is the 
optimal block size for this code. Another interesting class of 2D topological stabilizer codes are the 
color codes |10]. The color codes offer little practical advantage over the surface code if the goal is 
to build a quantum memory as some of the parity checks involve more than 4 qubits and the noise 
thresholds are similar as for the surface code, see e. g 



2.1 Surface Code 

The surface code is a version of Kitaev's toric code [1] in which the periodic boundaries of the torus 
have been replaced by open boundaries ^3] . Many of its properties and ideas for its use as a quantum 
memory were first analyzed in the seminal paper ^B]- The topological 2D realization of the CNOT 
gate was first proposed in |171 HE] . 

A simple sheet, depicted in Fig. [6} can encode one logical qubit. The linearly-independent parity 
checks are weight-4 plaquette Z-checks Bp and star X-checks As which mutually commute and are 
modified at the boundary to act on 3 qubits, see Fig. [6j Note that the star operators are just plaquette 
operators on the dual lattice when interchanging X ^ Z. The smallest surface code encoding 1 logical 
qubit which can correct 1 error is the code [[13,1,3]]Q Z is any Z-string which connects the north and 
south rough boundaries, we can deform this string by multiplication by the trivially-acting plaquette 
operators. X is any X string (on the dual lattice) connecting the smooth east and west boundary. 
As these strings have to connect boundaries in order to commute with the check operators, their 
minimum weight is L. Thus for general L the code parameters are [[L^ + {L — 1)^, 1, L]]. 

Using 13 qubits to correct 1 error does not seem very efficient, but the strength of the surface code 
is not fully expressed in its distance which only scales as the square root of the number of qubits in 
the block E3 

Let us consider how quantum error correction can take place first assuming that the parity check 
measurements are noise-free. If a single X error occurs on an edge in the bulk of the system, then 
the two plaquette operators next to it will have eigenvalue —1. The places where these plaquette 
eigenvalues are —1 are called defects. A string of X errors will thus produce two defects at its 
boundary. If the X error rate p is sufficiently small, one obtains a low density of close-by defects. Such 
errors are correctable as defects can be locally paired without much ambiguity. As we know, inferring 
an error E' which differs from the real error E by only stabilizer operators, plaquette operators in this 

^°One can minimize the qubit overhead while keeping the distance equal to 3 by rotating the lattice and chopping 
off some qubits to get a total of 9 qubits, see e.g. [49] . 

^^One can prove that the distance of any 2D stabilizer code is at most 0{L) 50 . However, one can also show [5T] 
that any block of size Rx R where R is less than some constant times the distance, is correctable, i.e. all errors in such 
Rx R patch can be corrected. These arguments show that there are no other 2D stabilizer codes with better distance 
scaling and that this scaling allows one to correct failed blocks of size beyond the distance. 
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Figure 6: Surface Code on a L x L lattice. On every edge of the black lattice there is a qubit, in 
total + (L — 1)^ qubits (depicted is L = 8). Two types of local parity checks, As and Bp, each act 
on four qubits (except at the boundary). The subspace of states which satisfy the parity checks is 
two-dimensional and hence represents a qubit. Z (yellow) is any Z-string connecting the north to the 
south boundary while X (blue) is any X-string connecting the east to west boundary of the lattice. 



case, is harmless. Here it means that we decode correctly as long as E'E does not represent X, an 
X-string which goes from one boundary to the other. From this picture it is clear that there should 
be a finite asymptotic threshold pc for noise-free error correction. 

For the bulk system the error syndrome thus gives us the location of the defects. A minimum- 
weight decoding algorithm then corresponds to finding a minimum-weight error string E{X) which has 
these defects as end-points. This decoding algorithm can be implemented using Edmond's minimum- 
weight matching (Blossom) algorithm [32] • Ideal decoding is not minimum weight decoding, but 



maximum-likelihood decoding as described in Section 1.1.1 One can estimate the threshold Pc as 



suming optimal maximum likelihood-decoding by mapping the problem onto a phase transition in a 
2D random-bond Ising model [53] . The value for pc assuming noise-free parity checks, ideal decoding 
and independent X errors is pc ~ 11% [IS]. 

This picture gets modified when the parity checks are inaccurate. A simple way to model noisy 
parity checks is to assign a probability q for the parity check outcome to be inaccurate while in 
between the parity checks qubits undergo X and Z errors with probability p as before. In practice, 
one would expect the parity check measurements to induce some correlated errors between the qubits 
of which we take the parity. For example, for the parity Z-check one may expect to cause additional 
qubit dephasing when more information than merely the parity is read out. 

As the parity check measurements are no longer reliable one needs to change their use as an error 
record. For example, a single isolated defect which appears for a few time-steps and then disappears 
for a long time is more likely due to a faulty parity measurement outcome than due to a string of X 
errors (which terminates at the isolated defect) that suddenly pops up and disappears without trace 
after a few time-steps. The strength of topological codes for sufficiently large L (as compared to using 
small codes and code concatenation) is that noisy parity checks can be dealt with by repeating their 
measurement as the additional noise which the parity checks produce on the code qubits is local and, 
at sufficiently low rate, correctable. 

Both minimum weight decoding and maximum likelihood decoding can be generalized to this set- 
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ting. We extend the lattice into the third (time) dimension |l6] , see Fig. [7j Vertical links, correspond- 
ing to parity check measurements, fail with probability q while horizontal links fail with probability 
p. In minimum weight decoding the goal is now to find a minimum weight error E which has vertical 
defect links, where the parity check is —1, as boundary, see Fig. [7j If the parity check measurements 
are ongoing, one needs to decide how long a record to keep before starting to infer the possible errors 
which depends on the failure probability q. In the simple case when q = p the record length is taken 
as L [53]. An analytical lower bound on the noise threshold for q < p is derived in [46] with the 
value Pc > 1.1%. Numerical studies in (using minimum weight-decoding) show a threshold of 
Pc ~ 2.9% for p = q. Below the noise threshold the logical error rate p(p, L) ~ exp(— k(p)L) where 
ti{p) ^ 0.8 — 0.9 at p = Pc/S [531 132] • AH these results have been obtained for toric codes, assuming 
periodic boundary conditions: one may expect results to be somewhat worse for surface codes [M] . 

The practical relevant question is how much overhead L is needed before one is in the scaling 
regime where the pseudo-threshold is close to the asymptotic threshold Pc{L) ^ pc? The pseudo- 
threshold for a small code such as [[13,1,3]] is very tiny, certainly no higher than 0.1% Even 
though small codes have worse performance than large codes they could still be used as testbeds for 
individual components and error scaling behavior |56j . 
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Figure 7: Picture from [JB]: ID cross-section of the lattice in space, and time. Grey links correspond to 
non-trivial —1 syndromes. Errors which could have caused such a syndrome are represented by black 
links. Horizontal black links are qubit errors while vertical black links are parity check measurement 
errors. Note that a possible error E has the same boundary as the grey defect links: a likely error 
E (in the bulk) can be found by looking for a minimum-weighted matching of the end-points of the 
grey links. 



Minimum-weight decoding with Edmonds' matching algorithm is a good method for decoding if 
our goal is to realize a quantum memory (with or without encoded Clifford group operations). As 
one never needs to physically do any correction (see the notion of Pauli frame discussed in Section 



1.4), the measurement record can be stored. This data record can then be processed at leisure and 
used to interpret a final or Mg- measurement on the qubits. The realization of such quantum 
memory will require that the record of parity check measurements is obtained at sufficiently high rate 



""^^Fig. 11 in 55j estimates a pseudo-threshold of 6.8 x 10 for the [[41, 1, 5]] surface code for an encoded transversal 
CNOT gate. This may be worse by a factor 2 from the pseudo-threshold for a single encoded qubit. 
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compared the error rate, since a low rate stroboscopic picture of the defects (even if they are obtained 
perfectly) could potentially miss the occurrence of a logical error. 

If our goal is to implement universal quantum computation, then the decoding algorithm will 
have to keep up with the logical error rate of the quantum computer as every time a non-Clifford 
gate is realized in the quantum circuit, one needs to know whether to do an additional logical Pauli. 
Researchers have developed faster renormalization-group decoders [571 EH] which process the defects 
using parallel processing over the 2D or 3D lattice in time O(logL). The idea of the simple decoder 
in ^58j which works for any D-dimensional stabilizer codes is to recursively match defects locally. For 
a 2D surface code with perfect parity check measurements one starts by dividing up the defect record 
into local clusters of 0(1) size. In each cluster the algorithm tries to find a local error which removes 
the defect. If a cluster contains a single defect for example, then no such local error can be found. 
Thus the next step is to enlarge the linear size of the cluster by a factor of 2 and reapply the same 
procedure on the leftover defect record. The decoder stops when no more defects are present or when 
one has a reached a certain maximum number of iterations r scaling as r = O(logL). For the toric 
code with perfect parity checks, Ref. [5H] has obtained a noise threshold of pc = 6.7% using this RG de- 
coder while the RG decoder in |57] achieves 9% (minimum-weight decoding via matching gives 10.3%). 



One can ask whether quantum error correction for the surface or other topological codes in D = 2 
or higher is possible by purely local means. Is there a non-zero noise threshold Pc when error correction 
has to be performed by a 2D noisy automaton which locally interacts with the surface code qubits 
|i6] . The dissipative correction procedure described in Section 1.3 is an example of a purely local 
error correction mechanism. Such local dissipative mechanisms will not be able to efficiently remove 
error strings in the surface code beyond the shortest length as an error string is only heralded by its 
two end-point defects. The local observation of one defect gives no information where the position 
of another matching defect is, except if one searches the entire neighborhood around the defect. One 
should contrast this with the 4D toric code where an error cluster is surrounded by defect surfaces 
such that the cluster can be removed by locally shrinking the defect surface. This weakness of the 
surface code is a generic property of 2D stabilizer codes and ties in with the discussion of self-correcting 
quantum memories which are unlikely to exist in 2D [50] . 

The classical non-local minimum- weight matching decoder can be viewed as a way of jump-starting 
our quantum memory with the power of reliably non-local classical computation. Note that the 
RG decoder is non-local (even allowing for parallel processing of clusters on the lattice by local 
automata) as the maximum number of recursions r scales as O(logL) leading to a maximum cluster 
size proportional to the linear size of the lattice. The lowest levels of the RG decoder are of course local 
and will provide some measure of protection by local means. In [21] the author devises an intricate 
scheme for doing purely local quantum error correction for the surface code by a 2D noisy cellular 
automaton assuming independent local errors; the mere existence of such a scheme is nontrivial. 

The idea of realizing the surface code in superconducting systems using ancilla qubits for measure- 
ment was considered in [59]. Ref. [56j consider experimental benchmarks for parity check measure- 
ments in a circuit QED setting. In order to reduce qubit overhead and possibly make better and more 
direct use of the cavity mode-transmon interaction, Ref. [60] considers a direct 3-qubit parity check 
measurement. In this scheme 3 transmon qubits are dispersively coupled to two transmission- line 
resonators such that the phase-shift picked up by an incoming microwave-pulse is (ideally) the same 
for qubit states with the same parity but differs by tt for qubit states with different parities. A simpler 
scheme for a two-qubit parity measurement in the cavity QED setting was proposed in Ref. [61] . 
In [62] the dispersive coupling of a set of transmon qubits to the cavity mode was proposed for the 
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measurement of a multi-qubit parity check. 



2.1.1 Multiple Qubits and Operations by Lattice Surgery 



How do we prepare the surface code memory in the states |0) , |l) or |±)? And how do we read out 
information, that is, realize and ? We prepare all elementary qubits in Fig. [6] in the state |0) 
and start measuring the parity checks with the goal of preparing the state |0). The state |00 . . . 0) 
has Bp = +1 and Z = +1 while the star operators Ag have random eigenvalues ±1 corresponding to 
the presence of many Z errors. Thus we choose some correction E for these Z errors (we pick a Pauli 
frame): the choice will not matter as E commutes with Z. If the preparation of |0) and the parity check 
measurements are noisy, one needs to measure the parity checks for a while before deciding on a Pauli 
frame for both X and Z errors. The preparation of |l) and |±) can be performed analogously using 
the ability to prepare the elementary qubits in |1) and |±) respectively. A projective measurement 
in, say, the |o) , |l)-basis (Mg-) proceeds essentially in reverse order. One measures all qubits in the 
Z-basis. Using the past record of parity Z-check measurements and this last measurement, one infers 
what X errors have taken place and corrects the outcome of Z = ±1 accordingly. 

There are various ways in which we can encode multiple qubits in the surface code and do a logical 
Hadamard and CNOT gate. The simplest and most efficient method in terms of qubit overhead is to 
encode multiple qubits in multiple separate sheets (as in Fig. |6]) laid out next to each other in a 2D 



array as in Fig. 10 Using operations on a single sheet one can do a logical Hadamard gate A 
CNOT between qubits in separate sheets can be realized using the idea of lattice surgery in which 
sheets are merged and split as proposed in |19]. One can also encode a qubit in an arbitrary state 
into the code, see [IHl |19] . 

We first consider the Hadamard gate for which HXH = Z and HZH = X. Consider doing a 
Hadamard rotation on every elementary qubit on a sheet. The resulting state is a +1 eigenstate of 
the Hadamard-transformed parity checks HAgW and HBpW which are the plaquette Z-check resp. 
the star X-check of the code Sduai defined on the dual lattice. The Hadamard gates map Z onto X^uai 
and X onto Zduai- In the dual code Sduai the rough and smooth boundaries are interchanged so that 
the lattice is effectively rotated by 90°. In order to complete the Hadamard gate and return back to 
the original code one can modify the parity checks at boundaries such that rough becomes smooth 
and smooth becomes rough again, see the details in |l6l HSj. 
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Figure 8: CNOT via 2-qubit quantum measurements. Here Mxx measures the operator X ® X etc. 
The ancilla qubit in the middle is discarded after the last measurement Mx disentangles it from the 
other two input qubits. Each measurement has equal probability for outcome ±1 and Pauli corrections 
(not shown, see Eq. ([T])) depending on these measurement outcomes are done on the output target 
qubit. 

The construction for the CNOT is based on the circuit in Fig. |8] which implements the CNOT gate 
through 2-qubit parity measurements, originally described in Ref. [63]. One could verify Figjsjby con- 
sidering the input states as a stabilizer subspace with stabilizer Z2 and logical operators Xi, Zi, X3, Z3, 
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using the Knill-Gottesman simulation technique, see Section 7 in [35] . One can also consider the evo- 
lution of the input \c)^ [0)2 1^)3 for bits c = 0, 1 and t = 0,1 exphcitly (here 1 denotes the top qubit 
in the Figure). For Mr^r^ = +1, we have a bit b^x = and Mxx = —1 corresponds to b^x = 1 etc. We 



have 



10)2 |t)3 ZtX', \c), ^(100)23 + 111)23) 

"^^^^ ZtXi \c),Xt \c), Xt |c)3 "^-^^ \c), l+)2 Ztxt \c © t)3 • (7) 

We observe the logic of the CNOT gate on qubits 1 and 3 in addition to corrective Pauli's Z^'^'^X^'"' 
which depend on the outcomes of the measurements Mxx and Mzz- The measurement Mx on the 
second qubit ensures that no information leaks to that qubit so that the CNOT gate properly works 
on any superposition of inputs. 

This circuit identity implies that we can realize a logical CNOT gate if we have the capability of 
projectively measuring the operators X ^X and Z ®Z oi two qubits encoded in different sheets. The 
capability to prepare a sheet in |0) and the measurement was discussed before. The realization of 
such joint measurement, say, X ® X ys, possible by temporarily merging the two sheets, realizing the 
measurement and then splitting the sheets as follows. Consider two sheets laid out as in Fig. [9] where 
a row of ancillary qubits is prepared in |0) between the sheets. We realize a rough merge between 
the sheets by including the parity checks, plaquette and star operators, at this boundary. If the 
parity check measurements are perfect, the new weight-4 plaquette Z-checks have +1 eigenvalue as 
the ancilla qubits are prepared in |0). The 4 new star boundary checks have random ±1 eigenvalues 
subject to the constraint that the product of these boundary checks equals the product of Xs of the 
two sheets. Hence a perfect measurement would let us do a X (g)X measurement. As the parity check 
measurements are imperfect, one needs to repeat the procedure in the usual way to reliably infer the 
sign of X ® X. We are however not yet done as we wish to realize a projective X ®X measurement 
on the qubits encoded in two separate sheets. This means that we should split the two sheets again: 
we can do this by reversing the merge operation and measure the ancillary qubits in the Z-basis and 
stop measuring the 4 boundary X-checks. Again, if the parity check measurements are perfect the 
plaquette Z-checks at the boundary of both sheets will be randomized, but both correlated with the 
outcome of the Z-measurement of the ancillary qubits. Hence the individual X eigenvalues of the 
separate sheets may be randomized, but they are correlated so that X ® X remains fixed. Similarly, 



a smooth merging and splitting (as between qubits C and INT in Fig. 10 ) with the ancillary qubits 



prepared and measured in the X-basis accomplishes a Z ® Z measurement. 



The procedure for a CNOT in Fig. 10 then consists of: first a preparation of the INT qubit in |o), 
then a rough merge and split of qubits T and INT followed by a smooth merge and split between 
qubits INT and C followed by a final Mj^ measurement of qubit INT. 



2.1.2 Topological Qubits and CNOT via Braiding 

A different way of encoding multiple qubits and realizing a CNOT gate was first proposed in HZ] . 
The disadvantage of this method is that it has an additional qubit overhead. A distance-3 smooth 
hole qubit (see the description below) costs many more than 13 physical qubits but such overhead or 
redundancy could also lead to greater robustness. 

In order to see how to encode multiple qubits we start with a simple square sheet with all smooth 
boundaries which encodes no qubits. Fig. [Il|^a). For a sheet with all smooth boundaries, the star 
operators are linearly-dependent while the plaquette operators are linearly independent so that no 
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Figure 9: Picture from [IH]: two sheets are merged at their rough boundary by placing a row of ancilla 
qubits in the |0) state at their boundary and measuring the parity checks of the entire sheet. For a 
similar smooth merge, the ancillary qubits in between the two sheets are prepared in the |+) state, 
see the INT and C sheets in Fig. 
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Figure 10: Picture from ^49j: using an ancilla (INT) qubit sheet we can do a CNOT between the 
control (C) and target sheet (T) by a sequence of mergings and splittings between the sheets. 



qubit is encoded. To encode qubits one makes a hole in the lattice, that is, one removes some checks 
from the stabilizer S. This is a change in topology which affects the code space dimension. In stabilizer 
terms: when we remove one plaquette, say, Bp^ for some p^, from the stabilizer S, then Bp^ is no longer 
an element in S but still commutes with S, therefore is a logical operator. The matching logical 
operator which anti-commutes with it starts at the hole and goes to the boundary. This encoded qubit 
has poor distance namely ci = 4 as i?p* is of weight 4. We can modify this procedure in two ways such 
that logical qubits have a large distance and its logical operators do not relate to the boundary. The 
particular choice of logical qubits will allow one to do a CNOT by moving holes. 

To get a logical qubit with large distance we simply make a bigger hole. We remove all, say, k'^ 
plaquette operators in a block (and all (fc — 1)^ star operators acting in the interior of this block) 
and modify the star operators at the boundary to be weight-3, no longer acting on the qubits in the 
interior, see FigjlT|(a). The qubits in the interior of the block are now decoupled from the code qubits. 
The procedure creates one qubit with Z equal to any Z-loop around the hole. The X operator is a 
X-string which starts at the boundary and ends at the hole. Clearly, the distance is the minimum of 
the perimeter of the hole and the distance to the boundary. We call this a smooth hole as the hole 
boundary is smooth. Of course, we could do an identical procedure to the star operators, removing a 
cluster of stars and smaller subset of plaquette operators and adapting the plaquette operators at the 
boundary. Such qubit will be called a rough hole and its X operator is a X-string around the hole (a 
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string on the dual lattice) and Z is a Z-string to the boundary. 




(b) 

Figure 11: (a) A smooth hole is created by removing a block of plaquette operators and the star 
operators acting on qubits in the interior of the block. The Z-loop around the hole is Z while X is 
an X-string to the boundary. The qubits inside the hole (4 in the picture) are decoupled from the 
lattice, (b) Two smooth holes can make one smooth qubit and two rough holes can make one rough 
qubit so that moving a smooth hole around a rough hole realizes a CNOT gate. 



In order to be independent of the boundary we use two smooth holes to define one smooth or 
primal qubit and use two rough holes to define one rough or dual qubit as follows. Consider two 
smooth holes 1,2 and define a new smooth qubit as |o)^ = |0,0)j^2 \^)s ~ |-'-'-'-)i2' -^^^ ^^^^ 
smooth qubit s we have Zs = Zi,i = 1,2 (we can deform Zi into Z2 by plaquette operators) and 



Xs = X1X2 which we can deform to an X-string which connects the two holes, see Fig. (ll)(b). The 
distance of this smooth qubit is the minimum of the distance between the holes and the perimeter of 
one of the holes (assuming the boundary is sufficiently far). Similarly, we can create a rough qubit by 
taking two rough holes and defining 

|0>. = ^ (|0'0>3,4 + 1 1' 1)3,4) ' |1>. = ^ (|0' 1>3.4 + 1 1' 0)3.4) • 

With this choice is the loop X3 (or equivalently X4) while Z^ = Z^Z^ is equivalent to the Z-string 
connecting the holes. 

Imagine moving smooth hole 1 around rough hole 3 as in Fig. [lT|(b). After the move, the X-string 
connecting the smooth holes will additionally go around the rough hole enacting the transformation 
Xg — )■ Xr ® Xs. Similarly, the Z-string Zr connecting the rough holes will wind around the smooth 
hole, leading to the transformation Zr — )■ Z s® Zr- The loops Z^ and X^ are not changed by the 
move. This action precisely corresponds to the action of a CNOT with smooth qubit as control and 
rough qubit as target 

The ability to do a CNOT with smooth qubit as control and rough qubit as target seems limited 
as all such gates commute. However one can use the one-bit teleportation circuits in Fig. |4]to convert 
a smooth qubit into a rough qubit and a rough qubit into a smooth qubit, using only CNOTs with 
smooth qubits as controls. We have already shown how to realize the other components in the one-bit 



^■^The action of the CNOT in the Heisenberg representation is Xc0lt — > Xc^Xt, Ic®Xt — > Ic<E)Xt, Zc®It 
and Zt ^ Zc® Zt where X^ {Xt) stands for Pauh X on control qubit c {t) etc. 
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teleportation circuit such as and M^. Thus by composing these circuits we can do a CNOT 
between smooth qubits alone (or rough qubits alone). 

How is the braiding done in practice? The advantage of realizing topological gates in stabilizer 
codes (as opposed to braiding of Majorana fermions or non-Abelian anyons in quantum Hall systems) 
is that braiding can be realized by changing where we measure the parity checks, or deforming the 
code. For example, one can enlarge the hole in Fig. 11 to include, say, 2 more plaquettes and 3 
more qubits in the interior. We stop measuring those two plaquette checks and the star checks in 
the interior, modify the star boundary measurements and measure the qubits in the interior in the 
X-basis. The modified weight-3 boundary checks will have random ±1 eigenvalues as their previous 
eigenstates were perfectly entangled with the qubits in the interior. This corresponds to a high Z-error 
rate around the modified boundary. By repeating the measurement and solidifying the results one 
can correct these Z-errors, but of course we may partially complete a Z-loop this way. The protection 
against a full Z-loop around the hole is thus provided by the part of the hole boundary which remains 
fixed. 

This implies that the hole can safely be moved and braided in the following caterpillar manner. 
One first enlarges the hole (while keeping its 'back-end' fixed providing the protection) so that it 
reaches its new position. In terms of parity check measurements it means that from one time-step to 
the next one, one switches from measuring the small hole to the large hole parity checks. Due this 
extension errors will occur along the path over which the hole is moved and if error correction is noisy 
we should not act immediately to infer the new Pauli frame, but repeat the new check measurements 
to make this new frame more robust. Then as a last step, we shrink the hole to its new position and 
corroborate the new measurement record by repetition. Alternatively, one can move the hole by a 
sequence of small translations, so that the hole never becomes large. The speed at which the hole 
can then be safely moved is determined by the time it takes to establish a new Pauli frame (eliminate 
errors) after a small move. Details of hole moving schemes are discussed in e.g. [HI], s 



2.1.3 Surface Code with Harmonic Oscillators 

It is possible to define a qudit stabilizer surface code, see e.g. [SS], where the elementary constituents 
on the edges of the lattice are qudits with internal dimension d and the code encodes one or several 
qudits. An interesting special case is when we take d ^ oo and each edge is represented by a harmonic 
oscillator with conjugate variables q. The goal of such continuous-variable surface code is to encode 
a non-local oscillator into an 2D array of oscillators such that the codestates are protected against 
local shifts in p and q. In addition, one can imagine using continuous-variable graph states to prepare 
such encoded states and observe anyonic statistics [66] . 

To get a surface code, we replace Pauli X by X{b) = exp(27ri6p) and Pauli Z by Z{a) = exp(27ria q) 
for real parameters a,b such that [Zi(a)Z^^(a), Xi(6)X2(6)] = for all a,b. In the bulk, a plaquette 
operator centered at site p can be chosen, for example, as Bp(a) = Zp^^io,) Z~^^{a) Zp^y(a) Z~^~{a) 
while Asib) = X;},{b)X,^^{b)X,_^{b)X;l^{b), see Fig. 
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Here Zp^x{,o) = exp(27rmgp_£) where qp-x is the position variable of the oscillator at site p — x 
{x and y are orthogonal unit- vectors on the lattice). One can observe from the Figure that Bp{a) 
and -Bp(a) commute with As{b) and Al{b) for all a,b in the bulk and at the boundary. We can 
thus define hermitian operators with real eigenvalues in [—1,1] as Hp{a) = l(_Bp(a) + B^{a)) = 
cos{27ia{qp-s:-qp+x + qp-y-qp+y)) and Hs{b) = ^{As{b)+Al{b)) = cos{27ib{-ps-x+Ps+x+Ps~y-Ps+y))- 
We now define the code space as the +1 eigenspace of all Hp{a), Hs{b) for all a and b. It follows that 
states in the code-space obey the following constraints: a state in the code-space is a delta-function in 
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ZXa) 

Z(a) Q Z"\a) X(b) 

Z(a) x^^^^ 

Z(c) 



X(b) 



X{d) 



Figure 12: Small example of the oscillator surface code where oscillators on the edges are locally 
coupled with plaquette and star operators so as to define an encoded oscillator with logical, non-local, 
displacements X{d) and Z{c). The realization of the four-oscillator interaction will require strong 
4-mode sqeezing in either position (at plaquettes) or momenta (at stars). 



the positions of the oscillators around all plaquettes p, that is, S^qp^x — Qp+x + Qp-y — Qp+y), while it is a 
delta-function in the momenta of the oscillators 5(— p^-x +Ps+£ +Ps~y ~Ps+y) located around all stars 
s. Perhaps one can compare such highly-entangled non-local state with the simpler cousin, the 2-mode 
EPR-state. In the 2-mode case we have two commuting operators Zi(a)Z^^(a) and Xi{b)Xi(b) on 
oscillator 1 and 2. The single state which is the eigenstate of cos(27ra(gi — ^2)) and cos{27ib{pi+p2)) 
for all a, b is the two-mode infinitely-squeezed state 6{pi + P2)S{(li — 12) ■ 

For this surface code, the operators Z{c) = exp(27rzc Ylie-yi where the path 71 runs straight from 
north to south commute with all Hp{a), Hs{b), see Fig. 12 Similarly, we have X{d) = exp{27rid Yljej2 
where 72 runs straight from east to west. As Z{c)X{d) = e"^^'^'^'^^^'^'^'^^ X (d) Z (c) , we can interpret Z{c) 
and X{d) as phase-space displacements of the encoded oscillator with logical position and momentum 
p = X]ie72 P'' ^ ~ Sie7i li- '^^^ deform these non-unique logical operators to follow deformed 
paths, e.g. multiply Z{c) by -Bp(c) plaquettes (note that if we multiply by Bp{c') with c' 7^ c we get 
an operator with the union of supports). 

How would one use such code to encode quantum information and what protection would it 
offer? As its qubit incarnation, a sufficiently-low density of independent errors on the lattice can be 
corrected. In reality however, each oscillator i will independently suffer from small shifts Zi{e)Xi{e') 
with |e|, |e'| <S 1. A local mechanism which enforces the perfect correlation constraint at plaquettes 
and stars can perform a local error correction, however the likelihood to make logical errors of the 
form Z{c)X{d) for small c,d will be high, as we are attempting to encode an analog signal. This 
of course is the drawback of any code encoding a continuous-variable. However, one can imagine 
using a 2-dimensional subspace of the encoded oscillator space, and use the qubit-into-oscillator 
code, see the end of Section L2, to prepare and use a discrete set of states in this subspace. In 
theory, one could prepare the encoded states |0) , |l) , |-|-) , |— ) of this surface code by preparing each 
local oscillator in the qubit-into-oscillator logical states |0) , |1) , |-|-) , |— ) and subsequently projecting 
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onto the perfectly correlated momenta and position subspace. For example, the state |0) of a local 
oscillator i is an eigenstate of Sg{a) = e^'^*'^'/", Sp = e~^*^''^ and the local Zi = This implies 

that after projecting onto the space with Hp{a) = 1, Hs{b) = 1 for all a, b, it will be an eigenstate of 
Z(l/(2a)) = e*'^^'67i ^'/"^ i.e. the encoded |o). The preparation of the codewords for this oscillator 
surface code would require considerable strength and flexibility of squeezing interactions. 



2.2 Bacon-Shor Code 



An interesting family of subsystem codes are the Bacon-Shor codes |67] which are competitive with the 
surface code for small numbers of qubits. For the [[n^, 1, n]] Bacon-Shor code the qubits are laid out in 
a 2D n X n square array, see Fig. [TSj The stabilizer parity checks are the double Z-column operators 
Z||^j for columns i = 1, . . . n — 1 and double X-row operators X=j for rows j = 1, . . . n — 1. It is also 
possible to work with asymmetric Bacon-Shor codes with qubits in a x m array. Asymmetric codes 
can have better performance when, say, Z-errors are more likely than X-errors (when T2 ^ Ti), see 



68] . The gauge group Q (see Section 1.1.1 ) is generated by weight-2 vertical XX links and horizontal 



ZZ links and contains the parity checks. The logical operators (which commute with Q but are not 
in S) are the single Z-column Z and a single X-row X. 




(a) 




(c) 



Figure 13: (a) 9-qubit [[9, 1, 3]] (Bacon) Shor code with black qubits on the vertices. Depicted are the 
2-qubit XX and ZZ gauge operators, (b) [[25, 1, 5]] Bacon-Shor code with X, a row of Xs (blue), and 
Z, a column of Zs (yellow). The stabilizer generators are double columns of Zs, Z|| j (one is depicted) 
and double rows of Xs, X= j.(c) Decoding for X errors (or Z errors in the other direction). Black 
dots denote the places where the double column parity checks Z|| j have eigenvalue —1 (defects). The 
X error string E has X errors in the red region and no errors elsewhere and Ec is its complement. 
Clearly the string E has lower weight than Ec and is chosen as the likely error. 
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(a) 



(b) 



Figure 14: (a) In order to measure the XX and ZZ operators one can place ancilla qubits (open 
dots) in between the data qubits. For example, for the ZZ measurement the ancilla qubit is placed 
in the |0) basis followed by 2 CNOTS on its adjacent data qubits to collect the syndrome and a final 
Z measurement, (b) Alternatively, to measure Z||^j one can prepare a 3-qubit entangled cat state 
;^(|000) + |111)) (green dots) which interacts locally with the adjacent system qubits. X=_i could be 
measured by preparing a cat state for the ancilla qubits at, say, the red dots. The ancilla qubits at 
the open dots can be used to prepare the cat states. 



Consider the correction of X errors sprinkled on the lattice, assuming for the moment that the 
parity check measurement of Z|| j is noise-free. For each column we note that an even number of X 
errors is a product of the vertical XX gauge operators and therefore does not affect the state of the 
logical qubit. This means that per column only the parity of the number of X errors is relevant. The 
double column operator Z|| j determines whether this parity flips from column i to column i + The 
interpretation of the eigenvalues of Z||^j is then the same as for a ID repetition code (or ID Ising 
model) with parity checks ZjZj+i. Double columns where Z|| j = ZjZj+i = —1 are defects marking 
the end-points of X-strings (domain walls in the ID Ising model). Minimum- weight decoding is very 
simple as it corresponds to choosing the minimum weight one between two possible X-error strings: 
E or the complement string E^. which both have the faulty double columns defects as end-points, see 



Fig. 13|^c). The code can thus correct all errors of weight at most [|J for odd n. Higher- weight errors 
can also be corrected as long they induce a low density of defects on the boundary. Note however that 
the number of syndrome bits scales as n whereas the number of errors scales with rt^. This means 
that in the limit n — )■ oo the noise-free pseudo-threshold Pc{n) — )■ as the fraction of incorrectable 
errors will grow with n. So, how do we choose n in order to minimize the logical error rate p{p,n): 
in [69] the authors find that the optimally-sized Bacon-Shor code for equal X and Z error rate p is 
given by n = ^ and for that optimal choice they can bound p{p,n) < exp(— 0.06/p). 

How do we obtain the non-local parity check values? One can either measure the XX and ZZ 
gauge operators and use this information to get the eigenvalues of X= j and Z|| j, or one measures 
the parity checks directly. The first method has the advantage of being fully local: the ancilla qubits 
for measuring XX and ZZ can be placed in between the data qubits, see FigjT4|a). In the second 
method we can prepare an n-qubit cat state, see e.g. [68]. We could measure Z\\^i using the circuit 
in Fig. [l](a) with a single ancilla in the |-|-) state and controlled-phase gates (CZ) gates. However, 
a single X error on the ancilla qubit can feed back to the code qubits and cause multiple Z errors 
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making the procedure non fault-tolerant. In addition, the interaction between the ancilla qubit and 
the code qubits is non-local. Instead, we encode the ancilla qubit |+) using the repetition code, i.e. 
we prepare the n-qubit cat state :^(|00 ... 0) + 1 11 ... 1)) such that a CZ-gate acts between one cat 
qubit and one code qubit. The ra-qubit cat state, which is stabilized by ZjZj+i and Xi . . . X„, can be 
made by preparing |+)®" and measuring ZjZj+i using local ancilla qubits. The ZiZi+i eigenvalues 
are recorded to provide the Pauli frame. In [68] further details of this scheme are given, including 
estimates of the noise threshold for asymmetric Bacon-Shor codes. 

Consider what happens when the local XX and ZZ checks are measured inaccurately. The good 
news is this only causes local errors on the system qubits. The bad news is that if the measurement 
outcome has some probability of error q, then the error probability for the nonlocal stabilizer checks 
outcomes will scale as qn. This is a disadvantage of the Bacon-Shor code. Previous work [TOl [68] has 
sought to improve the fault-tolerance of the parity check measurements by replacing the preparation of 
simple single qubit ancillas by fault-tolerant ones (methods by Steane and Knill). In [70] a best noise- 
threshold of Pc ~ 0.02% was numerically obtained for the (concatenated) [[25,1,5]] code. Recently, 
[n^ considered an alternative way of making the syndrome more robust, namely by simple repetition 
of the XX and ZZ measurements and a collective processing of the information (as is done for the 
surface code). We can view the effect of repetition as extending the ID line of defects to a 2D lattice 
of defects, as in Fig. [7] so that minimum weight decoding corresponds to finding a minimum weight 
matching of defect end-points. The error rate for vertical (black) links representing the parity check 
errors scales with n while the error rate for horizontal links (when one column has an even and the 
other column has an odd number of errors) scales, for low p, also with n. In |B2] the authors estimate 
that the optimal size for the Bacon-Shor code is then n ~ 0.014/p and that for this choice, the failure 
probability p{p) < exp(— 0.0068/p). Hence for an error rate of p = 5 x 10^'^ we can choose n = 28 
giving a logical error rate of p ~ 1.25 x 10^^. 



2.3 Subsystem Surface Code 

One can ask whether it is possible to reduce the requirement to measure 4-qubit parities by using sub- 



system codes ^ The Bacon-Shor code exemplifies this, but the Bacon-Shor code is not a topological 
subsystem codes as the stabilizer checks are nonlocal on the 2D lattice. Several topological subsystem 
codes have been proposed in which weight-2 parity checks are measured [72], but the asymptotic noise 
threshold for such codes is typically quite a bit lower than for the surface code, see e.g. [73] . 

In [H] a topological subsystem code was proposed, — a subsystem surface code — , in which the 



non-commuting parity checks are of weight-3 and the stabilizer generators are of weight 6, see Fig. 15 
More precisely, the gauge group Q is generated by the triangle operators XXX and ZZZ, including 
cut-off weight-2 operators at the boundary. The stabilizer group S = QnC{Q) is generated by weight-6 
plaquette operators (at the boundary — )■ weight-2 operators). By measuring, say, the Z-triangles we 
can deduce the eigenvalues of the Z-plaquettes which are used to do error correction. 

For a. L X L lattice one has a total of 3L^ + 4L + 1 qubits and 2L^ + 4L independent stabilizer 
generators which gives + 1 qubits. One of these qubits is the logical qubit whose Z and X commute 
with all Z and X-triangles. Similar as in the surface code, a vertical Z-line through 2L qubits can 
realize Z while a horizontal X-line realizes X. The logical operators for the gauge qubits, one 
for each plaquette, are pairs of triangle operators on a plaquette generating the group Q. One can 
multiply, say, the vertical Z-line by Z-triangles to become a Z which acts only on L qubits: in 



""^^One can prove that 2D qubit codes defined as eigenspaces of 3-local mutually commuting terms are trivial as 
quantum codes |7T] . 
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it is indeed proved that the distance of the code is L. Note that such weight-L Z acts on the logical 
qubit and the irrelevant gauge qubits. 

For a code with distance L = 3 one thus needs 41 elementary qubits, substantially more than 
for the surface code. Multiple qubits can be encoded in this subsystem code by making holes as 
for the surface code. One can expect that braiding and lattice surgery methods for this code can 
be established in the same way as for the surface code. The interesting feature of this code are its 
relatively-high noise threshold obtained by reduced-weight parity checks (at the price of a bit more 
overhead). Decoding of stabilizer syndrome information is done by interpreting the syndrome as 
defects on a virtual lattice which can processed, similar as for the surface code, by minimum weight 
matching of defects or by RG decoding. For noise-free perfect error correction and independent X, Z 
noise, the authors report a maximum threshold of pc ~ 7% (compare with 11% for the surface code). 
For noisy error correction the threshold depends on how single errors with probability p in the parity 
check circuit affect the error rate on the virtual lattice. Modeling this effective noise-rate on the 
virtual lattice, the authors find a noise threshold of Pc ~ 0.6%. 

It is not surprising that decoding for the subsystem surface codes can be done using the decoding 
method for the surface code. It was proved in [75] that any 2D topological subsystem or stabilizer 
code can be locally mapped onto copies of the toric code. The upshot is that for any such code one 
can find, after removing some local errors by local correction, a virtual lattice with toric code parity 
checks and an underlying effective error model. An example of another 2D topological subsystem 
code which may be analyzed this way is a concatenation of the [[4, 2, 2]] code with the surface code. 
If we use the [[4, 2, 2]] code as subsystem code then the concatenated code has weight-2 and weight-8 
check operators. The scheme may be of interest if the weight-2 checks can be measured fast and with 
high accuracy. 

A possible lay-out of the elementary qubits for the subsystem surface code with superconducting 
transmon qubits coupled through resonators (including ancilla transmon qubits for measurement) was 
considered in |76j. The idea of a direct 3-qubit parity measurement [60] via a parity-dependent shift 
of the cavity mode in the circuit QED setting could be useful in reducing the need for ancillas, in 
addition to being making better use of the given physical interactions. The authors of [71| report a 
threshold at least 0.97% by modeling such direct parity measurement. 

3 Discussion 

Physics in the past century has demonstrated the experimental viability of macroscopic quantum 
states such as the superconducting state or Bose-Einstein condensates. Quantum error correction 
which strives to preserve not a single macroscopic quantum state but the macroscopic states in a 
small subspace can be viewed as a natural but challenging extension to this. At the same time storing 
macroscopic quantum information is a first step towards the more ambitious goal of manipulating 
quantum information for computation purposes. 

The current qubit realizations seem perhaps awkwardly suited to constitute the elementary qubits 
of an error-correcting code. Most elementary qubits are realized as nondegenerate eigenlevels (in a 
higher- dimensional space), approximately described by some Hq = ^Z. The presence of Hq immedi- 
ately gives a handle on this qubit, i.e. processes which exchange energy with this qubit will drive it 
from |1) to |0) and vice versa (Rabi oscillations) and coupling of the qubit to other quantum degrees 
of freedom can be used for qubit read-out. Passive (non-time dependent) interactions with other 
quantum systems are intentionally weak and only lead to significant multiple-qubit interactions if we 
supply energy in the form of time-dependent AC fields meeting resonance conditions. To drive, keep 
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Figure 15: Picture from [71]: subsystem surface code on a lattice of size L x L with square 
plaquettes (depicted is L = 3). The qubits hve on the edges and vertices of the plaquettes and are 
acted upon by weight-3 X and Z-triangle operators (which are modified to become weight-2 operators 
at the boundary). The stabihzer checks are weight 6 except at the boundary. 

or project multiple qubits via local parity checks in a code space where they are highly entangled, 
active control at the elementary qubit level will thus be continuously needed, making the macroscopic 
coding overhead look daunting. 

A different route toward the realization of quantum error-correction is based on passive Hamilto- 
nian engineering. Instead of obtaining the desired stabilizer interactions dynamically through tem- 
porarily coupling qubits, one may attempt to engineer a physical system such that it has efffective, 
say, 4-body Hamiltonian interactions of the surface code between nearby qubits in a 2D array [77] . 
The difficulty with this approach is that it requires one, two and three-body terms in the effective 
Hamiltonian to be small: the elementary qubits should therefore have approximately degenerate levels 
|0) and But, then in order to encode information in the ground-space of such Hamiltonian, one 
will need to hft this degeneracy to address these levels. A possible advantage of this approach is that 
the protection is built into the hardware instead of being imposed dynamically. For example, the 
dynamical simulation of a topologically-protected Hamiltonian in a rotating frame where the strong 
single-qubit terms are absent does not guarantee that the physical system will also be driven to the 
ground-space of this simulated Hamiltonian; this will depend on (the time-scale of) the coupling with 
the dissipative environment. In contrast, topological protection at the hardware level means a realiza- 
tion of a Hamiltonian together with its low-temperature physical environment. A disadvantage of this 
passive approach is that the desired 4-body interactions will be typically be arrived at perturbatively 
and thus their strength and therefore the gap of the topologically-ordered Hamiltonian versus temper- 
ature may be small. An example is a 2D array of superconducting rhombi in a uniform magnetic field, 
each rhombus consisting of four superconducting islands, which can provide a Hamiltonian realization 
of the Bacon-Shor code with XX check and non-local stabilizer Z|| j terms; see the review [75] for 
this and various other protected-qubit proposals for superconducting qubits. Another example is the 
proposal to realize the parity checks of the surface code through Majorana fermion tunneling between 
2D arrays of superconducting islands, each supporting 4 Majorana bound states with fixed parity 
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[79] . One drawback of this Hamiltonian-based 'self-correcting' approach is that a finite temperature 
'quantum memory phase' may only exist in 3D or 4D quantum systems, see e.g. the review [80]. 
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